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1 Introduction 

Suppose that we are selecting n points, pi,P2, ■ ■ • ,Pn, at random in a rectangle, say R = [0, 1] x [0, 1] (see 
Figure nl). We denote by tt the configuration of n random points. With probability 1, no two points have 
same x-coordinates nor y-coordinates. An up/right path of tt is a collection of points pi-^,pi^, ■ ■ ■ , pi^. such that 
x(Pii) < x(pi2) < ■ ■ ■ < x(pi^) and yipii) < y(,Pi2) < ■ ■ ■ < y(Pik)- The length of such a path is defined by the 
number of the points in the path. Now we denote by ^„(7r) the length of the longest up/right path of a random 
points configuration tt. 
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Figure 1: random points in a rectangle 

As one can see from Figure y, a configuration of n points gives rise a permutation. For the example at 
hand, the corresponding permutation is (51324)- Therefore we can identity random points in R and random 
permutations, and we use the same notation tt. In this identification, /„(7r) is the length of the longest increasing 
subsequence of a random permutation. 

The longest increasing subsequence has been of great interest for a long time (see e.g. |l|, p9], M). Especially 
as n ^ 00, it is known that E(;„) - 2^ ||l, |9|, ^ (also % H, H) and Var(;„) ~ con 
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some numerical constant cq = 0.8132 • • • . Moreover, the limiting distribution of /„ after proper scaling is 
obtained in Q in terms of the solution to the Painleve II equation (see Section ^ for precise statement). 
An interesting feature is that the above limiting distribution function is also the limiting distribution of the 
(scaled) largest eigenvalue of a random GUE matrix |^ , the so-called "GUE Tracy- Widom distribution" F2 . 
In other words, properly centered and scaled, the length of the longest increasing subsequence of a random 
permutation behaves statistically for large n like the largest eigenvalue of a random GUE matrix. There have 
been many papers concerning the relations on combinatorics and random matrix theory : we refer the reader 
to H, |I|, H, |l|, H, § I |6|, 0, H H |3|, |, |3|, |2^, |, §. The purpose of this paper is to survey the analytic 
results of the recent papers [g| ^ and discuss related topics. 

In random matrix theory, three ensembles play important roles, GUE, GOE and GSE (see e.g. Pql). Since 
random permutation is related to GUE, it would be interesting to ask which object in combinatorics is related 
to GOE and GSE. For this purpose, we consider symmetrized permutations. In terms of random points, 5 
symmetry types of the rectangle R are considered, denoted by the symbols D, 0, S, □, and S. Throughout 
this paper (and also in p, 0]), the symbol ® is used to denote an arbitrary choice of the above five possibilities. 
Let S ~ {{t, t) : < t < 1}, the diagonal line, and (5* = {{t, 1 — t) : < i < 1}, the anti-diagonal line. Consider 
the following random points selections : 

(i). n : select n points in R at random. 

(ii). : select n points in R\d and m points in d at random, and add their reflection images about S. 

(iii). S : select n points in R\d* and m points in S* at random, and add their reflection images about 5*. 

(iv). : select n points at random in R, and add their rotational images about the center (1/2, 1/2). 

(v). : select n points in R\ 6, m+ points in S and ?7i_ points in 6* at random, and add their reflection 
images about both S and S*. 

Define the map t on S'„ by i{x) = n + 1 ~ x. Let fp(7r) denote the number of points satisfying n^x) — x, and 
fpi(7r) denote the number of points satisfying 7r(a;) = l{x) (fpi represents negated points : see Remark ^ below). 
Each of the above process corresponds to picking a random permutation from each of the following ensembles 

^n^m = {tT e S2n+m ■ TT = 7r"\ fp(7r) = Ul} (1.2) 

Sn,m = {tT e S2n+m ■ TT = ITT^^t, fpi(7r) = m} (1.3) 

S^ ^ {n e S2n ■■ TT ^ lttl} (1.4) 

Sn,r,i+,rn- = {tT G 54„+2m++2m_ : TT = TT^S TT = tTT^^t, fp(7r) = 2to+, fpi(7r) = 2m_}. (1.5) 

We denote the length of the longest increasing subsequence (equivalently, the longest up/right path) of tt in 
each of the ensemble respectively by 



Remark 1. The map tt i-^- i^^Tr gives a bijection between 5*^^ and S^^. Thus L^ „, has the same statistics with 
the length of the longest decreasing subsequence of a random involution with ra fixed points taken from S^^. 
From the definition, L^ „ is the random variable describing the length of the longest increasing subsequence of 
a random involution taken from the same ensemble. 

Remark 2. We may identify S2n with the set of bijcctions from {—n, • • • , —2, — 1, 1, 2, • • • ,71} onto itself. In this 
identification, S^ becomes the set of signed permutations ; 7r(a;) = — 7r(— x). The longest increasing subsequence 
problem of a random signed permutation is considered in [p6| and ||7| . 

Remark 3. Under the identification in Remark El S^^ ^_ becomes the set of signed involutions with m+ fixed 
points and to_ negated points (we call x a negated point if ■k{x) ~ —x.) 

In this paper, we are interested in the statistics of L® as n —^ 00. Especially for 0, S and Kl, we are 
interested in the cases when m = [\/2na] for [Z\, m ~ [V2n/3] for S, and to+ ~ [y/na] and m^ = [\/nf3] for 
S with fixed a, f3 > where [k] denotes the largest integer less than or equal to k. Then for most cases, the 
expected values have the same asymptotics. Namely, if we set N = n,2n + m,2n + m, 2n, An + 2m+ + 2to_ for 
each of D, 0, S, Q, Kl case respectively, we have 

li,n 5^ . 2, (L7) 



when < a < 1 and /3 > are fixed for 0, S and Kl. When a> 1, we have different expected value in the limit 
(see Section ||.) 

On the other hand, the variance behaves asymptotically like coiV^/'^ but now with different constant cq 
depending on the symmetry type. It is because each symmetry type has different limiting distribution : L^ has 
GUE fluctuation, L^„^ GOE fluctuation and L^ GUE^ fluctuation (see Section ^ below for precise statements). 
Here GUE^ denotes the statistics of a superimposition of eigenvalues of two random GUE matrices. Similarly 
for GOE^. The cases of and show more interesting features. For 0, the limiting distribution function 
changes depending on the value of a = r7i/v2n. The fluctuation is GSE when a < 1, GOE when a = 1 and 
Gaussian when a > 1. By taking suitable scaling limit a ^ 1, we can find a certain smooth transition between 
GSE and GOE. For 0, the value a = m^/ ^/n determines the limiting distribution ; the value 7n_ plays no role 
in the transition. The fluctuation is GUE when a < 1, GOE^ when a ~ \, and Gaussian when a > 1. 

In Section y, we define the Tracy- Widom distributions for GUE, GOE and GSE as well as new classes of 
distribution functions describing the transition around a = 1. Main results are stated in Section 0, and Section 
4 includes some applications and the related problems. Most of the results in this article are taken from [g, ||. 



The only new result is Theorem 4.2 
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2 Tracy- Widom distribution functions 

Let u(x) be the solution of the Painleve II (PII) equation, 



with the boundary condition 



Uxx — 2w + xu, 



u(x) ^ —Ai(x) as X — > +cx), 



(2.1) 



(2.2) 



where Ai is the Airy function. The proof of the (global) existence and the uniqueness of the solution was first 
established in 11^ : the asymptotics as x ^ — oo are (see e.g. |l^, O]) 



u{x) = -Ai{x)+0 



^{4/3)2 



,1/4 



as x ^ +00, 



u{x) 



1 + 



(;?))• 



as X 



Recall that Ai{x) 



2y^a;i/4 



as X 



-oo. Define 
v{x) := 



{u{s))'^ds, 



(2.3) 
(2.4) 

(2.5) 



so that v'{x) = {u{x))'^. 

We introduce the Tracy- Widom distributions. (Note that q := —u, which Tracy and Widom used in their 
papers, solves the same differential equation with the boundary condition q{x) ^ +Ai(x) as a; — > oo.) 

Definition 1 (Tracy- Widom distribution functions). Set 



F{x) := expl - / v{s)ds 
E(x) := expl - / u{s)ds 



expl-i 



(s — x){u{s)Y ds 



and set 



^2(X) 

F^{x) 



Fix) 



F{x)E{x) 



exp[ — / {s — x){u{s)) ds], 

{F^{x)Y'\^-^r<^)'^^, 



Fi{x) := F{x) [E{x)-^ + E{x)] /2 = (^2(2;)) 



1/2 



3-^ ir «(s)ds ^ gi /~ u{s)ds 



/2. 



(2.6) 

(2.7) 

(2.8) 

(2.9) 
(2.10) 



In |3^ and [Q, Tracy and Widom proved that under proper centering and scaling, the distribution of the 
largest eigenvalue of a random GUE/GOE/GSE matrix converges to -^2(2:) / Fi{x) j F^(x) as the size of the 



matrix becomes large. We note that from the asymptotics ( p.3[) and (2^), for some positive constant c, 

^3/2 N 



F{x) = l + 0(e- 

E{x) = l + 0(e""^' 

F(x) = 0(e-"l^l') 

E(x) = 0(e-=l-l''') 



as a; ^ +00, 
as a; — > -l-oo, 
as a; ^ —00, 
as a; — > —00. 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



Hence in particular, HtjIx^+oo Ffj{x) = 1 and HtjIx^^oo Ffsix) — 0, /3 ~ 1,2,4. Monotonicity oi Fp{x) follows 
from the fact that Fp{x) is the limit of a sequence of distribution functions. Therefore Fi3{x) is indeed a 
distribution function. 

As indicated in Introduction, we need new classes of distribution functions to describe the phase transitions 
from xgse to xgoe and from xgue to Xgoe^- First we consider the Riemann-Hilbert problem (RHP) for the 
Painleve H equation [Q Q. Let T be the real line M, oriented from +oo to — oo. Let m{- ; x) be the solution 
of the following RHP : 

m{z; x) is analytic in z G C \ F, 

1 _g-2i(|z^+x2;) 

to(z; x) = I + O(-) as z — > cxd. 
Here m^{z]x) (resp,. rri-) is the limit of m{z';x) as z' -^ z from the left (resp., right) of the contour F 



m-^-{z; x) = TO_(z; x) 



for z e F, 



(2.15) 



m±{z]x) — \mi^iam[z =p ic,x). Relation (2.15) corresponds to the RHP for the PII equation with the special 
monodromy data p = —q = l,r = (sec H, EOi, also M, M). In particular if the solution is expanded at 

Z — CXD, 



we have 



m(z; x) = / + ^ ' + 0{-^ 



2i(rni(a;))i2 = ~2i(mi(a;))2i 
2i(rni(a;))22 = -2i{mi(x))ii 



as z ^ oo, 
= u{x), 



(2.16) 

(2.17) 
(2.18) 



where u{x) and v[x) are defined in (2.1)-( p.5| ). Therefore the Tracy- Widom distributions above are expressed 
in terms of the residue at oo of the solution to the RHP ( p.l5| ). It is noteworthy that the new distributions 
below which interpolate the Tracy- Widom distributions require additional information of the solution of RHP. 



Definition 2. Let m{z]x) be the solution of RHP (2.15) and denote by mjk{z\x) the (j/i;)-entry of m{z;x). 
For w > 0, define 

F^(x;w) :— F{x)< [m22{~'iw; x) — ■mi2{—iw]x)]E{x)^^ + [m22i—iw;x) + mi2{—i'w;x)]E{x) >/2, (2.19) 



(2.20) 



and for w < 0, define 

F^{x;w) ■.^ei'"'-^'""F(x) 

X < [— m2i(— iw; x) + mii(— iw; a;)]£'(a;)~"'" — \m2i{—iw\x) -\- mii{—iw]x)\E{x)> /2. 
Also define 

m22i~iw;x)F2{x), w > 0, 



F^ix;w) 
F^ix;w) 



-el^'-'^^'"m2i{-iw;x)F2ix), w < 0. 



(2.21) 

(2.22) 



First F^{x;w) and F^{x;w) are real from Lemma 2.1 (i) below. Note that F^{x;w) and F^{x;w) are 
continuous at w = since at z = 0, the jump condition of the RHP ( p. 15 ) implies (mi2)+(0; x) = — (mii)_(0; x) 
and (m22)+(0;x) = -(m2i)-(0; a:). In fact, F^{x;w) and F^{x;w) are entire in w e C from the RHP ( |2.15D . 

From (|2.1l[) - (^.14|) and (^.24|) -( ^.27|) below, we see that 



lim F^{x;w),F'^{x;w):^l, lim F^(a;; w), F^(x; u;) = 



X — >■ — oo 



(2.23) 



for any fixed w € M. Also Theorem 3.3 below shows that F^{x;w) and F^{x;'w) are limits of distribution 
functions, implying that they are monotone in x. Therefore, F'^{x;w) and F^{x;w) are indeed distribution 
functions for each w e M. 

We close this section summarizing some properties of m{~iw; x) in the following lemma. In particular the 
lemma implies that F'^{x;w) interpolates between F4{x) and Fi{x), and F^{x;w) interpolates between -^2(2^) 
and Fi{x)'^ (see Corollary O). 

Lemma 2.1. Let 0-3 = ( i*! ) , cri = ( ? J ) , and set [a, b] ~ ab — ba. 

(i). For real w, m{~iw; x) is real. 

(ii). For fixed w G K., we have 



m{—iw; x) 



')) 



'1 _p^w^ — 2xw^ 



m{-zw;x) = (^+(e-^"'''))(_^_a!3+2.. J 



w > 0, X ^ +00, 



w < 0, X ^ +00, 



(2.24) 
(2.25) 



m{—iw; x) 
(Hi). For any x, we have 



m{-ivu;x) ^ 1 ^-iV(-|-^+--k3g(-^(-x)^/^+y2»^(-x)^/^V3^ w>0,x-.-oo, (2.26) 



V2 



-1 1 



w <0, x^ -00. (2.27) 



lim mi—iw.x) — lim aim(—iw,x)ai = , 



10^0+ 



'E{xf' 
E{xf ^ 



(2.28) 



(iv). For fixed w G R\ {0}, m{—iw]x) solves the differential equation 

d 



dx 



m = w[m, a^] + u{x)(7im, 



(2.29) 



where u{x) is the solution of the PII equation (2J.), ( |2.2| ) 



Corollary 2.2. We have 



(2.30) 
(2.31) 
(2.32) 
(2.33) 
(2.34) 
lim i^^(x;w) = 0. (2.35) 



F^{x;0) = 


= Fi{x), 


lim F {x;w) = 


-- Fi{x), 


lim F {x;w) = 

— > — oo 


-- 0, 


F^(x;0) . 


- F,{xf 


lim F {x;w) - 


-- F2{x), 



in — > — oo 



Proof. The values at w = follow from ( 2.28 ). For w —> ±oo, note that from the RHP ( 2.15| ), we have 



Imiz^oo'iT^iz'jx) = L n 

3 Main Results 

As in the Introduction, let N denote n, 2n + m,2n + to, 2n, An + 2m+ + 2to_ for each of D, 0, S, □, Kl case 
respectively. We scale the random variables : for permutations and involutions, 

□ _ Lg - 2^/]V _ Ll„, ~ 2Vn ^ _ Lg^„, ~ 2Vn 

and for signed permutations and signed involutions, 

1^-2^ Ll^,,m^ - 2Vn 

All the results in this section are taken from [g| which utilizes the algebraic work of [Q . 

First, we state the results for random permutations and random signed permutations. Random permutations 
show GUE fluctuation in the limit, while random signed permutations have GUE^ fluctuation. 

Theorem 3.1. For fixed x e M, 

lim Pr(x°<a;)=^2(x), (3.3) 

n — >OQ 

lim Pr(xg <x) ^F2{xf. (3.4) 

n — *oo 

For the involution cases, we have the following limits. 
Theorem 3.2. For each fixed a and (3, and for fixed x G R, we have : for 0, 

„l™,PK<[v^„] < ^) = ^4(0.), < a < 1, (3.5) 

„l^P^«[v^l^-)=^i(^)' (3-6) 



forS, 

i™ P^X^^I^,] < -) = ^i(-), /? > 0, (3.8) 
and for K, 

^lim Pr(x,^[^„j f^^] <x)^ F^ix), < a < l,/3 > 0, (3.9) 

Jim Pr(x,t[^„],[^^] < ^) = 0, a > l,/3 > 0. (3.11) 

This theorem shows that for and K, the Umiting distributions differ depending on the value of a. As 

indicated earher in the Introduction, as a ^ 1 at a certain rate, we obtain smooth transitions. From Corollary 



2.2, the following results are consistent with Theorem 3.2 



Theorem 3.3. For fixed w eR, f3 > and x eR, 

Mm Pr(x^^<x) ^F^{x;w), m = [V2^ - 2w{2n)^/% (3.12) 

lim Pr(x^„^ m_ < a:) - F^{x;w), m+ = [^-2wn^^^], m_ = [V^/3]. (3.13) 

When a > 1, Theorem |3.2| shows that we have used inappropriate scaling. In a proper scaling, we obtain 
normal distribution A^(0, 1). 

Theorem 3.4. For fixed a > 1 and f3 > 0, as n ^ oo, 



L0,/_ .-{a + l/a)VN 
L^ ~(a + l/a)VN 



A^(0, 1) in distribution, (3.14) 



A^(0, 1) in distribution. (3.15) 



v/(l/a-l/a3)iVi/4 

All the above results are on the convergence in distribution. We also have convergence of moments for all 
the cases. 

Theorem 3.5. For each case of the above theorems, all the moments of the random variable converge to the 
moments of the corresponding limiting distribution. 



From this result, we can obtain the asymptotics of variances. Especially for D, the variance is 

hm^ 



hm ^^^^^ = r x^dF2ix) - ( / xdF^ix) ) . (3.16) 



oo 



Evaluating the integrals, we obtain the value 0.8132 • • • (see \p7\). 

The outline of the proofs is as follows. First we consider the Poisson generating function. It is to let the 
number of points be Poisson. For example, we define 

^.°(i) = e-*^^i^Pr(L°<0. (3.17) 

n=0 



The de-Poissonization lemma [E4J tells us that in the limit n ^ oo, Pr(L^ < l) ^ -F/°("-^); hence it is enough 
to obtain the asyniptotics of the Poisson generating function. The crucial point is that there is a determinantal 
formula for each Poisson generating function. For the cases D and □, the determinant is of Toeplitz type 1 18, 31|, 



while for the rest, it is of Hankel type Q. In fact, as in Q, there are general identities between sum of Schur 
functions and determinantal formulae (see |5| for details) , which can be used to consider other type of Young 



tableaux problems (see Subsection 4.4 below). Now general theory connects Toeplitz/Hankel determinants and 
orthogonal polynomials. It turns out that to analyze all the above cases, only one set of orthogonal polynomials 
are needed, namely the orthogonal polynomials on the unit circle with respect to the weight e^* "^"^ ^ d9 / {2tt) . 



Now following Fokas, Its and Kitaev 1 15|, there is a Riemann-Hilbert representation for orthogonal polynomi- 
als. Let E be the unit circle in the complex plane oriented counterclockwise. Let Y{z) be a 2 x 2 matrix-valued 
function satisfying 



zGE, (3.18) 




as z 



where Y+{z) (resp. Y-{z)) denotes the limit of Y{z') as z' ^ z satisfying \z'\ < 1 (resp. \z'\ > 1). Then 
one finds that for example, the 11 entry of Y{z) is the fc"' monic orthogonal polynomial with respect to the 
weight e'^*'^°^^dO/{2'K). Once we have a Riemann-Hilbert representation, we can employ the steepest-descent 
method (Deift-Zhou method) developed by Deift and Zhou ||l3| to find asymptotics as parameters become large 
(or small). For our case, the parameters are t and fc, and taking proper scaling, we obtain precise asymptotics 
which eventually yield the convergence in distribution and convergence of moments. We also mention that in 
our analysis, equilibrium measures play a crucial role as in the papers [|ll| y, |lO| . 

4 Applications and related topics 

4.1 Random involutions and random signed involutions 

The ensemble S^^ is the set of involutions with n 2-cycles and m 1-cycles. In the previous section, we considered 
the limiting statistics when n and ra are related by to = ['\/2na] with a being finite ; either fixed or a — *■ 1 with 
certain rate. It is of interest to consider the whole set of involutions without constraints on the number of fixed 
points. Similarly, the signed involutions without constraint on the number of fixed points and negated points 
is also of interest. We define the ensembles of involutions and signed involutions 

^„ = {7re5„:7r = 7r-i}, (4.1) 

S*^ = {tT e S-an : TT = TT-l, TT = tTTi}, (4.2) 

and denote by L„ (tt) and i^ (tt) the length of the longest increasing subsequence of tt G Sn and that of tt e S^ , 
respectively. 



Noting Sn = U2fc+m=n Sf,m^ ^6 have 

Pr(Z„ < ;) = 



1 



Sk 



J2 Pr{Ll„^<l)\S, 



k.m I 



(4.3) 



2fc+Tn— n 

It is not difficult to check that (see pp. 66-67 of [E6|) as ti ^ oo, the main contribution to the sum |5fc_,„| 



I 'S'fe m I comes from 



/2fc - (2fc)"+i/'' < m < V2fc + (2fc)"+i/''. Comparing with the scaUng 



Z_-/2fc-t-Tn— n 

[v2fc — 2w{2k)^''^] in Theorem 3.3, the main contribution to the sum (4.3) comes from when w = 0, or a = I. 



Thus we obtain GOE Tracy- Widom distribution function in the hmit. Similarly, we can obtain the convergence 
of moments. The signed involution case is analogous. 



Theorem 4.1. For fixed x E 



lim Pr 

n — *C30 



2y/E 



n 



1/6 



< X 



( L^-2^ 

^'SL^'y 22/3^1/6 <^ 



^Fi{x), 
= F,{xf. 



(4.4) 
(4.5) 



We also have convergence of all the moments. 

This result should be compared with the results on random permutation and random signed permutation 
where the limiting distribution was F2{x) and F2{x)'^ under the same scaling of the above (see y], [Q, and 



Theorem 3.1 above). 

Remark. As in the permutation case, the length of the longest increasing subsequence and the length of the 
longest decreasing subsequence of random involutions have the same statistics. This can be seen by noting that 
there is a bijection (called Robinson-Schensted correspondence, see e.g. [gg) between the set Sn of involutions 
of n letters and the set of standard Young tableaux of size n, and the rows and the columns of standard Young 
tableaux have the same statistics under the push forward of the uniform probability distribution on Sn under 
this bijection. 



4.2 /9-Plancherel measure on the set of Young diagrams 

Let Yn be the set of Young diagrams, or equivalently partitions, of size n. Given a partition A = (Ai, A2, • • • ) h n, 
let d\ denote the number of standard Young tableaux of shape A. We introduce the /?-Plancherel measure M^ 
on Yn defined by 



mP{\) 






AGr„. 



(4.6) 



When j3 — 2, this is the Plancherel measure which arises in the representation theory. We are interest in the 
typical shape and the fluctuation of A where A is taken randomly from the probability space F„ with M!^. 

A motivation introducing the above measure is the result of Regev [Q . In [Q , it is proved that for fixed 
/3 > and fixed I, as n ^ 00, 






{d.y 



f/2lr. 



(V27r)('-i)/2n('-i)('+2)/4 






i/3ii:-?-Q 



Xk 



\pd\ 



(4.7) 



i<k 



10 



The multiple integral on the right hand side is the Selberg integral which can be computed exactly for each (3. 
When (3 — 1,2, 4, this integral is the normalization constant of the probability density of eigenvalues in GOE, 
GUE and GSE, respectively (see e.g. |2^). So the basic question is if the f5 in the definition of the /3-Plancherel 
measure corresponds to the f3 in the random matrix theory. 

The well-known Robinson-Schensted correspondence [Q establishes a bijection between Sn and the pairs 
of standard Young tableaux with the same shape of size n, RS : tt i— > (P(7r), Q(7r)). Especially, we obtain 
Suhn'^M ~ 1*^"! ~ ""■■ Moreover, under RS, the length of the longest increasing subsequence of tt S 5„ is 
equal to the number of boxes in the first row of -P(7r) (or equally of Q{tt)). Therefore under RS, the Plancherel 
measure M^ is simply the push forward of the uniform probability measure on S'„ to Yn, and the number of 
boxes in the first row of a random Young diagram and the length of the longest increasing subsequence of a 
random permutation have the same statistics : GUE ((3 ~ 2) fluctuation in the limit. 

If RS(7r) = (P, Q), then RS(7r~-^) = (Q, P) (see e.g. Q). Therefore the set of involutions Sn is bijective to 
the set of (single) standard Young tableaux, and the number of boxes in the first row of a random Young diagram 
taken under the probability M]^ has the same statistics with the length of the longest increasing subsequence 
of a random involution : GOE (/3 = 1) fluctuation in the limit. 

In fact, it is shown in [^0[ || |l| that, for the case of /3 = 2, in the large n limit, the number of boxes in the 
fc*'' of a random Young diagram corresponds to the fc*'' largest eigenvalue of a random GUE matrix. Also the 
typical shape of A is obtained in |2^ and [g9) which is related to Wigner's semicircle law [g^. On the other 
hand, the second row of a random Young diagram for /3 = 1 is discussed in [pj implying that it corresponds to 
the second eigenvalue of a random GOE matrix. It would be interesting to obtain similar results for general 
row for /3 = 1 and also for general f3. 

4.3 Random turn vicious walker model 

Random permutations and random involutions arise also in certain random walk process in a one dimensional 
integer lattice. We call a particle left-movable (resp. right-movable) if its left (resp. right) site is vacant. 
Initially, p particles are located at the points 1,2,3, ••• ,p. At each time step t (t = 1,2,---), one particle 
among the left-movable particles is selected at random, and is moved to its left site (so at i = 1, the leftmost 
particle is moved.) Suppose this process is repeated for n time steps. It is found that p7| there is a bijection 
between all the possible configurations and the set of Young tableaux with at most p rows of size n : simply, in 
the k^^ row, write the times when the particle is originally at the point k (see Figure g). In this correspondence 
the number of moves made by the k*"^ particle (counted from the left) is the number of boxes in the fc*^ row of 
the corresponding Young tableau. If we take the limit p —i- oo, we remove the constraint on the number of rows. 
In other words, we start with countably many particles located at 1, 2, 3, • • • and at each time step, we move a 
randomly selected left-movable particle to its left site. Then the statistics of the number of moves made by the 
leftmost particle in n time steps is identical to that of the number of boxes in the first row of a random Young 
diagram under the probability M^ (or the length of the longest increasing subsequence of a random involution) . 
Hence in the large n limit, we obtain GOE fluctuation. 

Now suppose after n steps of left moves, by taking n steps of right moves, we want the particles brought 
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Figure 2: random turn vicious walkers and Young tableaux 

back to their original positions. The first n steps gives a Young tableaux and the next n step gives another 
Young tableaux. So we are in the situation of pairs of Young tableaux, i.e. /3 = 2 case. Especially if the number 
of particles are infinite, the statistics of the moves made by the leftmost particle is identical to the length of 
the longest increasing subsequence of a random permutation : GUE fluctuation [ll^ . 



4.4 Symmetrized versions of Johansson's model 

In p2|, Johansson introduced a certain model which has several probabilistic interpretations, as a randomly 
growing Young diagram, a totally asymmetric one dimensional exclusion process, a certain zero-temperature 
directed polymer in a random environment or as a kind of first-passage site percolation model. Here we consider 
the symmetrized versions. The limiting distributions arc parallel to the results of Section ^ depending on the 
symmetry type. 

Each model is defined on M, a subset of 7?, and at each site {i,j) G M, we define a random variable w{i,i) 
of geometric distribution. We are interested in "the length of the longest up/right path". Let g{q) denote the 
geometric distribution with parameter q. Let < g < 1, and let a,/? > such that ay/q,(3y/q < 1. Otherwise 
stated, the random variables w{i,j) are independent each other. We denote by (j, fc) /^ {j',k') the set of 
up/right paths vr from (j, fc) to (j',k'). 



(i). Let M = Z2_. Let 



Define 



wiij) --g(<?). 



G°(7V)=max{ ^ u;(z, j) : ^ € (1, 1) / (TV, TV)}. 



(4.8) 



(4.9) 
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(ii). Let M = 1?. Let w{i,j) = if i = or j = 0, and for i,j ^ 0, let 

w{i,j) ==^w{-i,~j) ^ g{q). (4.10) 

Define 

G^{N)^mSix{ Y^ w{i,j) -.n e{-N,-N) /" {N,N)}. (4.11) 

{i,j)e-iT 



(iii). Let M = Z^,. Let 



w{hj) ^w{j,i) -- g{q), i^j, (4.12) 

w{i,i) -- 9{a^/q)- (4.13) 



Define 



G^(7V)=max{^ ii;(z, j) : ^ e (1, 1) / (^, ^)}. (4.14) 

(ij)e7r 

(iv). Let M = Z+ X Z_. Let 

w(i,-j) ==«;(.?, -i) - g{q), ij^-j, (4.15) 

^«(^,-^) ^ gWV^)- (4.16) 

Define 

G^(7V)=max{^ «;(*, j) : ^^ G (1, -A^) / (A^, -1)}. (4.17) 

{i,j)eTT 

(v). Let M = Z^. Let w(i, j) = if i = or j = 0. Otherwise 

w{i,j)=w{-i,-j) ^ g{q), \t\^\j\, (4.18) 

w{i,i) ^w{-i,-i) - g{a^), (4.19) 

u;(^,-^)=«;(-^,^) ^ .gC/^V?)- (4-20) 



Define 



G^(7V) =max{ ^ u;(i, j) : tt G (-iV, -TV) / (iV,iV)}. (4.21) 

(i,j)e7r 



We are interested in the fimiting distributions of G® (N) . 
Theorem 4.2. Set 



'yC'?) = -I ^ , p(<?) = — ^ — ^ — • (4.22) 
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I ■ 



(ii). 



(Hi). 



.--(^^^l^^#^-)-^(^)- 



( G^iN;q)-rjiq)i2N) 



lin, Pr('^!fe4-^Mi^<x) = f4W, 0<a<l, (4.25) 



W^oo \ p{q)N'^/^ 



G^{N;q,a)~T^{q)N ^ ^\ ^0,.,.., „ , 2w 



^^'^J^ piq)mir' ^-} - P^i-M. «-l-^(^)]^> (4.26) 

^^M^^^^^4.rj^<^) - 0, a>l. (4.27) 



w^oo \ p{q)N^/^ 



(iv). 



(v). 



Mm p,iG-iN;q,P)-,iq)N ^ ^^^^^^ ^ ^ ^^ ^^^^.^ 



Af^oo \ p{q)N^/^ 



f G^{N;q,a,f3)^rjiq)i2N) 
nZo'^^ 22/3p(9)(27V)l/3 

f G^{N;q,a,P)-r,{q)i2N) 
N^^'^^ 22/3p(g)(27V)l/3 

/ G^(iV;g,a,/3)-r;(g)(2JV) 

A^"?L"^ 22/3p(9)(27V)l/3 



i^2(2;), 


u;), 


o< 

a = 
a > 


a<l,/3>0, 
2w 


,/3>0 


(4.29) 
(4.30) 
(4.31) 


r [X, 
0, 


^ p(g)(27V)i/3 
l,/?>0. 



Remark. The results ( [4.2S| ) and ( 4.26 ) with w ^ are obtained in |£l[. For D, the longest up/right path in 
(1, 1) / (M, N), M ^ N, is also considered. 

The proof of the above theorem is analogous to that of the result in Section ^. There again is a Toeplitz/Hankel 
determinantal formula for each case. In fact, the only change with the symmetrized permutation case is that 
we have {1 + q + 2 ^ cos 9)^ dO / {2n) instead of e^*' ^°^ ^ dO / {2tt) for the weight. In analyzing the asymptotics of 
orthogonal polynomials, we need different scaling, but once we scale, the analysis is parallel. One important 
common property of the above two weights is that both of the supports of their equilibrium measures change 
from the full circle to a part of the circle with one gap depending on 2t/k < 1 and 2t/k > 1, and N/k < ri{q)^^ 
and N/k > ri{q)^^. And also in the one gap case, the equilibrium measures decay like a square root at the ends 
of the supports. 
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